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Abstract 

The energy dependence of such characteristics as a ratio of the total number of charged particles to the 
total flux of EAS Cherenkov radiation, a ratio of E'^/j^ > 1 GeV muon flux density at the distance of 600 m 
from a shower core to charged particle flux density, a ratio of the energy transferred to the electromagnetic 
component of EAS to the primary particle energy is presented. Their comparison with two-component mass 
composition of cosmic rays (p-Fe) in the framework of calculations by a QGSJET model is given. 

I. INTRODUCTION tion of breaks is in the detailed study of different 

EAS characteristics in the region of the first and 



second irregularities in the spectrum. We have 
made such a work in [|8l|9l[l0llIIl[l2l|Il. 



II. EAS CHARACTERISTICS IN THE SU- 
PERHIGH ENERGY REGION 



The irregularities in the cosmic ray (CR) en- 
ergy spectrum of "knee" type at E'o — 3 x 
10^^ eV and "ankle" type at ^ 8 x 10^^ eV 
found in [dl |2l are yet of special interest from 
the point of view of interpretation of these phe- 
nomena from the position of astrophysics. In 

recent years a few papers H 13 have been Longitudinal development 
published which try to explain such a behavior 

of CR spectrum with the help of new models The cascade curves of EAS development 
of generation and propagation of CRs. There in Figjl] were reconstructed according to the 
exists also another namely a nuclear-physical method suggested in [|l4l|. It is seen from Figjl] 
point of view for the formation of irregularity that the maximum depth of cascade curves de- 
at -Eq — 3 X 10^^ eV [[61 13. In some sense, the pends on the primary particle mass composi- 
answer to the problem on reasons of the forma- tion as well as the hadron interaction model. It 

is seen from Figjl] that to describe the experi- 
mental cascade curve (Xmax, Nq) the QGSJET 
model is better- suited [fTSll . So we use this 



* s.p.knurenko@ikfia.ysn.ru 


1 tema @ ikfia.ysn.ru 






^ i.ye. sleptsov.ikfia.ysn.ru 





10' 



10 













QGSJET(Fe); 


SIBYLL (Fe) X 




DPMJET(Fe); 


QGSJET (p) 




SIBYLL (p); 


DPMJET (p) 




• Yakutsk 





f I I 

200 400 600 800 1000 1200 1400 1600 1800 



X , [ g / cm' 1 

Figure 1: Comparison of the experimental cascade 
curve of HAS development (Eq = 10^^ eV) with 
different models of hadron interactions. 

model for the estimation of mass composition 
of primary particles. Fig j2] presents the calcu- 
lations of the maximum depth Xmax using the 
QGSJET model for the primary proton and iron 
nucleus, and experimental data obtained at the 
Yakutsk EAS array. It is seen that the velocity 
of shift of Xmax to sca level depends on the en- 
ergy range. In the framework of the QGSJET 
model the experimental data are indicative of 
the change of mass composition of primary par- 
ticles in the energy range £'o = 3xlO^^ — 3x 
10^*^ eV and at > 3 x 10^*^ eV. 

Radial development 

Fig J2] and Figj3]present experimental data: a) 
the density of muons with Eth^ > 1 GeV at a 
distance of 1000 m from the shower core [[TT|. 
b) the root-mean- square radius Rm.s. of charged 
particle (LDF) [fT2ll . These data are compared 
with calculations using the QGSJET model. 
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Figure 2: p^(lOOO) vs Eq relations for observed 
events 10^^ - 10^^ eV (squares), 10^^ - lO^o eV 
(points) and 10^° — 10^^ eV (triangles). Expected 
iblcr bounds for the distributions are indicated for p, 
Fe and 7 primary by different curve as in the legend. 

Protons, iron nuclei and 7-quanta are consid- 
ered as primary particles. A confidence interval 
taken in calculations is ±lcr. From Fig. 3 it fol- 
lows that in the energy range 10^^ — 5 x 10^*^ eV 
within the boundaries of confidence there are 
~ 60 % of showers for the iron nucleus and 
~ 90 % for the proton. At > 5 x 10^^ eV 
the portion of iron nuclei decreases, the per cent 
of protons and 7-quanta increases. It is seen 
from Fig j3] that the "heaviest" mass composition 
is observed at Eq ~ 10^^ eV. 

Correlation of EAS parameters 

The most sensitive instrument for the model 
of hadron interactions and mass composition of 
primary particles is the EAS muon component. 
At the Yakutsk complex EAS array the muons 
with Ethr > 1 GeV are measured by a shower 
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Figure 3: Dependence of the root- mean-square LDF 
radius of charged particles on energy. The curves are 
calculations by the QGSJET model for the primary 
proton and iron nuclei. 



registration in 70 % cases. Figj4] presents the 
correlation of Nfj, — Ng parameters and Fig|5] 
gives the portion of muons (ratio of muons to 
all charged particles) at a distance of 300 m 
and 600 m from the shower core. In the same 
place the calculations by the QGSJET model 
are given. The tendency for an increase of 
light nuclei in the cosmic ray primary flux at 
Eq > 3 X 10^^ eV is marked by these data. 

Energetic EAS characteristics 

At the Yakutsk complex EAS array the 
shower energy is determined by measurements 
of the total flux of EAS Cherenkov Ught F, 
the total number of charged particles, Ns, and 
muons, A^^, at sea level 11161 . Fig j6] presents 
the energy-dependence of Ng/F ratio. It is seen 
from calculations that the ratio strongly depends 
on a mass composition. From comparison of the 
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Figure 4: Total number of charged particles A^^ and 
muons Nf^i at sea level. The curves are a calculation 
by the QGSJET model for the primary proton and 
iron nuclei. 




Figure 5: Portion of the muons with Ethr > 1 GeV 
(%). p^(300)M(300) and p^(600)/p,(600). 



calculations by the QGSJET model for the pro- 
ton, iron nucleus and experimental data it fol- 
lows that the mass composition changes just af- 
ter the first knee in the spectrum, i.e. in the in- 
terval i?o = 5 X 10^^ — 10^'^ eV. It is evident from 
Fig jV] where the portion of energy transmitted in 
to the electromagnetic EAS cascade is shown. 



io°r™T- 



10" 



Yakutsk (small array) 
-QGSJET (A=1) 
QGSJET (A=56) 



QGSJET01, p(70%), Fe(30%) 
Yakutsk 8x10"< E <2x10"eV 
<E>=1*10"eV, n=857 
QGSJET01 p(100%) 
CNO(100%) 
Fe(100%) 



10' = 



10'° 10" 
E„ [eV] 



10'" 




Figure 6: The dependence of Ng/F on a primary 
energy. 



Figure 8: Fluctuations of the parameter Xmax at 

Eq ~ 10^*^ eV. 
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Figure 7: A portion of the energy transferred to the 
electromagnetic FAS component by Cherenkov light 
data at the Yakutsk array. 

Fluctuations of some EAS parameters 

In this section we consider the fluctuations 
of Xmax and Rm.s. obtained by the measure- 
ment of Cherenkov EAS light and density of 
charged particle flux (see Figj8] and Figj9]). In 
both cases dX Eq = 10^® eV fluctuations are 
considerable and correspond to the mixed mass 
composition of primary particles. From com- 
parison with calculations by the QGSJET model 



Figure 9: Fluctuations of Rm.s. at Eq = 10^^ eV. 
Curves are the QGSJET model calculations for the 
primary proton and iron nucleus. 

we have the following relationship: the light nu- 
clei are ~ 70 % and heavy nuclei are ~ 30 % 
(see the dotted line in FigJS]). 

III. RESULTS AND DISCUSSION 

At the Yakutsk complex array for almost the 
35-year period of continuous observations the 
unique experimental data on electron, muon and 
Cherenkov EAS components in the region of 
superhigh and ultrahigh energies have been ac- 



cumulated. Results on longitudinal and radial 
development of EAS are presented in figures. 
From the total combination of data one can se- 
lect two energy regions where, as is seen from 
Figures, the characteristics of EAS have the 
complex dependence on the energy. This is the 
energy ranges 1 - 1 eV and 1 - 1 eV. 
As is known, in these energy ranges the irreg- 
ularities of "knee" and "ankle" in the energy 
spectrum of EAS are observed. 

From the comparison of all experimental 
data with calculations by the QGSJET model 
in ifTTl the results of the cosmic ray mass com- 
position in the energy range 10^^ — 3 x 10^^ eV 



have been obtained. It follows from the anal- 
ysis that after the "knee" the mass composi- 
tion becomes heavier and in the region of "an- 
kle", on the contrary, becomes lighter. Such 
a conclusion doesn't contradict the hypothesis 
on a cosmic ray generation up to the energy 
~ 10^^ eV in our Galaxy and their propaga- 
tion according to the model of anomalous dif- 
fusion in the fractal interstellar medium. Begin- 
ning with Eq > 3 X 10^^ eV the mass compo- 
sition becomes lighter and it doesn't contradict 
the presence of cosmic rays metagalactic origin 
in the total flux of cosmic rays. 
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